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ABSTRACT 



This thesis presents some general equations describing the 
particle density distributions v»ithin a plasma. Solutions to 
these equations are discussed for various geometries, with particular 
emphasis given to the effects of diffusion losses at probes of dimen- 
sions comparable to or greater than the mean free path and to magnetic 
field end effects. For mapy of the geometries discussed, including 
one probe geometry, the effects of a superimposed magnetic field is 



discussed 



I 

INTRODUCTION 

One of the complex problems in the field of plasma physics is 
that of predicting the effects of surfaces vith vtiich the plasma is in 
contact upon the particle density distribution functions and related 
phenomena. No plasma is infinite, and consequently a plasma must be 
confined by vralls, which may be either conductors or insulators. 

Further, in order to make measurements within a plasma, probes may 
be inserted. Vhat is their effect on the particle density distri- 
bution and related phenomena such as potential distribution and 
energy dissipation? bhat are the particle currents to various surfaces? 
The latter are particularly important because they determine the 
heat liberated at the surface, which in turn is an energy loss which 
may cause very high temperatures, increased local ionization, and 
melting of materials in extreme cases. This thesis discusses these 
effects when the surfaces are insulators. 

In this thesis some general equations describing a plasma are 
presented vith their corresponding restrictions. Solutions to 
these equations are discussed for various geometries, with particular 
emphasis given to the effects of probes and to magnetic field end 
effects on the steady-state density distributions and on the energy 
losses at the surfaces involved. For many of the geometries dis- 
cussed, including one probe geometry, the effects of a superimposed 
magnetic field is studied. 
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II 



GENERAL RELATIONS 

In this section, general equations describing a plasma are 
developed and discussed. 

A. Conservation Equations 

The Boltzmann Transport Equation 

Knx)-Vr'(^w) -n('2:-V^)X'n(C)fcX v)*\5x =J}(x'-x')F(v)4<vy4^jeV<^ V 



is used as the starting point for the development of the macroscopic 
relations governing the plasma. By ascribing appropriate values to X 
and >• £13 [2], ve obtain the continuity equation; 



1^ +V-(ioYj^ = n0,.-vj ( 2 - 2 ) 

(in which y- is the mean ionizing frequency and the mean recombination 
frequency) , and the equation for conservation or momentum; 



in which the right hand term is really a definition of the mean 
collision frequency, >£ > the pressure *p is a symmetric tensor 

(2»3a) 

and in vrhich 




(2«3b) 



and 

( 2-3c) 

E, B, q, m, and vj are, respectively, the electric field intensity, 
the magnetic field intensity, particle charge, particle mass, and 
drift velocity. 
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The continuity equation ( 2 - 2 ) can be further simplified. V/e 
shall consider a tvo component system of electrons and singly-charged 
ions in tjhich volume ionization by electrons is essential. Recombina- 
tion takes place predominately at the v^alls, vdth volume recombination 
being negligible. V'e have, then, 

vhere V is the net ionizing frequency and the particle current density 
is 



i ^ ( 2— l:a) 

In the above equations, and in the remainder of this thesis, the 
subscript "+" refers to singly-charged positive ions, and the sub- 
script refers to electrons. Under those conditions for vtiich 
the rates of decay of ions and electrons are essentially equal, so 
that 

St ^ ^ 5t. 

t-T€ obtain for the continuity equation 

v-r. 

In the momentum equation (2-3), the pressure term is generally 
a difficult one to evaluate because *p is not known, and the 
dependence on its divergence makes good approximations difficult. 
However, for the simple case in vhich the distribution of random 
velocities is essentially isotropic, so that 



(2-5) 

(2-6) 



• < 







I 






m 













I 








the pressure teiisoi' reduces to 






n m v ' 



=->7bT = 



YieV 

M > 



( 2 - 8 ) 



where D and ^ are the diffusion and nobility coefficients respectively. 
The last tern is derived from the preceeding ones by the use of 
Einstein’s relations 



JI _ kr *2^'' 
~ 2e 



(2-9) 



This generalisation (equation (2-8)) is valid only for a plasma which 
is characterized by a small drift velocity, 

( 2 - 10 ) 

The last tern in equation (2-3) is frequently written 

V, := _ £ £ (2-11) 

which defines the mobility, yu . In the case of ion-neutral particle 

interactions, is a constant because the forces are inverse fifth 

power polarization forces. In all other cases v»e must introduce an 

appropriate average value, V , which we will consider independent 

of V ,. 
d 

One important consequence of the abov'e assumptions is that the 
ion and electron nobilities and diffusion coefficients are scalar 
constants. 

Substitution of equations (2-8) and (2-11) into the momentum 
equation (2-3), and rearranging terms, yields, for the positive 
ions of mass Ms 

xS •= - f], (2-12) 

and for the electrons (mass m); 

i yu_/l>^S = “ fl. 

h 



(2-13) 



For stationary (steady-state) situations in the plasma, vre find 



that these equations reduce to 

Q + X 0) — 



(2-lii) 



B. Current Density, Field and Potential Equations 
From Poisson’s equation vre bsve 



y-E = i: 



<Ca J 



(2-15) 



and from Maxvrcll's equations for the stationary caset 

VxE 



( 2— l6) 




These three equations, with the continuity equation (2-6) and the 
simplified momentum conservation equation (2-lii) e^qsress the 
governing conditions for a plasma vhich complies with the assumptions 
ennumerated to date. These will be solved for current equations and 
for general field and potential relations. 

From equation (2-6) it is readily deduced that 



= r. -t-(vx f) (2-i8) 

in vhich (V x F ) is a vector vhich we shall nov' determine 
Let us define three orthogonal directions in the plasma, one parallel 
to the magnetic field, designated by the subscript '* || ", one 
perpendicular to the field and to any surfaces parallel to the field, 
and the third transverse to the field and to the perpendicular 
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direction "T" 



In cylindrical coordinates, vrith a field parallel to 



the z-axis and cylindrical boundaries parallel to the field these 
vould be the z~, r^ and directions respectively. The ion and 
electron particle currents at the sheath edge are equal at all points 
for a steady-state plasma because they must arrive in equal numbers 
to recombine. For a geometry in vhich the insulator surfaces coicide 
with coordinate surfaces (the only geometries for vhich ve will be 
able to obtain exact solutions for the particle density distributions). 



This requirement is met either hy a plasma vhich has a rotational 
axis of symmetry parallel to the magnetic field or one in vhich the 
field lines coincide vdth one set of orthogonal coordinate lines and 
in vhich the resulting distribution is a function of only one coordinate. 
The latter condition vould be met in the vicinity of a plane surface 
parallel to the m.agnetic field lines in an "infinite" plasma. Here, 
for a y-z plane and a field parallel to the z-axis, the perpendicular 
direction v^ould be the x-direction and the transverse direction the 
y-direction. If these symmetry conditions are imposed, it is readily 
deduced that equations (2-19) apply throughout the plasma if the 
vialls are insulators. If, however, the end surfaces are conductors 
(i.e,, electrodes), equation (2-19b) does not apply. The parallel 




(2-19a) 



and 




(2-19b) 



at the sheath edges. We novr impose the condition that 



V-T 'A+ ^ ^ - Et- = O. 



( 2 - 20 ) 



6 



I 







a 



I 



r 



I 



I 

I 

I 



current is then a conduction current, not a diffusion current, and is 
not treated by the developments to follovT. This current does not 
affect the perpendicular currents, hov^ever, and such a plasma vill 
be treated as an infinite plasma in the longitudinal direction, vith 

( 2 - 21 ) 

From equation (2-lli) and the symmetry conditions (2-20), it is 
immediately seen that the transverse component equation iss 

rt--p = (2-22) 

Froiri this equation, and equations (2-19)> we can see that ^ ^ ) 
of equation (2-l8) has only a transverse component: 

(vxF) 

The parallel and perpendicular component equations derived from 
equation (2-lli) are 



and 



Vi\\ = 






(2-2lj) 



(2-25) 



A simple solution is possible if the assumption of proportionality 
holds: 

Vn,. ^ 

n, ~ n- n • 

Ey eliminating from equation (2-2h), v»e obtain 



( 2 - 26 ) 



n, = -DiV|,n_ - -Dm,.9;in- 



( 2 - 27 ) 



vtiere Dn|| is a quasi -ambipolar diffusion coefficient for diffusion 



7 



parallel to a magnetic field and Ds is the ambipolar diffusion 
coefficient which is commonly developed for a plasma when B ® 0 by 
infering congruence, 



r. -r.-r 



( 2 - 26 ) 



from the current density divergence equality equation (2-6) and from 
the proportionality condition (2-26) [23 Dl • Dg is related to 
Schotti^*s ambipolar diffusion coefficient, 



Do. - 






by the equations 



D< 






( 2 - 29 ) 



( 2 - 30 ) 



- D.(i •<- ^). 



(2-30a) 



In equation (2-30a) we have employed the following relations for 
charge density and conductivity 

p ^(Yl+-yi-)e (2-31a) 

and 

O' - ^ /U-yl^ , ^ 2-3 lb) 

substituting equation (2-22) into equation (2-25), ve obtain 



rV j_ + t Hti. = - D± , ( 2-32) 



By the use of equation (2-19a), may be eliminated from the 
electron equation, to yield 




I 






1 * 




1 I 




I 




( 2»3li) 



By simplifying and solving for Q , we obtain 



n = 



- DsNZlTJ- 



I -f 



I 4- 



0“ 



B"’ 






wnere Dmj^ is a quasi -ambipolar diffusion coefficient for diffusion 
perpendicular to a magnetic field. Thus diffusion perpendicular 
to the magnetic field is reduced by an inverse dependence on the 
square of the magnetic field intensity. 

If there is no superimposed electric field, so that E ® Es, 
v4iere Es is the space-charge field, ESj_ nay be obtained by substitu cing 
equation (2-3ii) into (2-32); 

_ 0 T), 



h- 



y\ 



(2-3?) 



Similarly, 



F- - P- " ^11 d _ ~ P- Mm (2-36) 

^1' /U— VI - n 

The transvei'se space charge field is zero as a consequence of our 

assumptions that there is no tr-ansverse density gradient. 

For a quasi-neutia 1 plasma, in wiiich 

(2“37) 

it follovjs from equation (2-30) that Ds * Da, and that 



E 



> 



D^- D- 

/tt+ + 



Vi.n 
n . 



(2-38) 



In the absence of a magnetic field, equation (2-3?) xtien simplified 
and added to equation (2-36), yields the single equation for the 
space-charge field 



Ds-D- % 
/A- n 



^ D.-q 

yM, +/<- VT ' 
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(2-39) 



From the field equations, we obtain the follovring equation for 

\ 

the space charge potential for a plasma vrith only a perpendicular 
density gradients 




In the absence of a magnetic field, equation ( 2~39) yields 



(2-liO) 




D.-Ds 







D- - D-*- 

^4 ir 




(2-U) 



The zero for potential in both of these equations is taken as the 
"center" of the plasma; that point (or region) where the density 
gradient is zero. 



C. The Diffusion Tensor 

The current density equations (2-22), (2-27) and (2-3ii) can be 
combined into one equation 



ri = 

vhich, in rectangular coordinates, is 



(2^2) 




(2-h3) 



and in cylindrical coordinates is 



( 
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(2-Ui) 



The components of the diffusion tensor are 



Drt|i = l>s, 


( 2“liSa 


"D/n - - 


(2-lt5b) 


LV f r| 1 — • j V 




(2-li5c) 


rVi^^T ^ ^ B> lXy\^ ^ 


(2-U^d) 



We shall now examine each of these components. Equation 
( 2 “U 5 a) reflects the fact that the magnetic field does not affect 
diffusion in the direction parallel to itself. Further, we have 
shown that 

'^^'’,1 - (2^1i6a) 

C't >1. ^ 

Allis ^23 has also shovm that 

Qs (2~U6h) 

^co 

where ^ is the conductivity at the center of the plasma. The quasi- 
neutral limit ( 3 “i 46 a) is applicable to a plasma in which the 
mobilities (and consequently the conductivity) are so great that 
negligible space charge can develop. The above limits are the lo\f 
and high limits for the ambipolar diffusion coefficient. 

The perpendicular component is more complex. First, as the 
magnetic field approaches zero, 






( 2-lt7a) 
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1 




As B becomes very large, DtHj^ approaches zero as 

P5 



B-yc^o 

From equation (2-Ii6a) it is seen that 






i2-hlh) 



D, 






IX 



■p->0 , 



(2-li7c) 



and from (2-60a) that 



Do,, 



D_g^ 



(2.1i7d) 



"1 0^ 6 3 ^ ^ 

viiich can be non-zero only if 0, As shorn in equation (2-^5) 

for a long, narrow plasma this can be true only if 



D 



m, 



■D_ 



X (3^-vcJ / 

From these equations it is apparent that for the parallel and 



(2-l47e) 



perpendicular diffusion tensor components in a magnetic field the high 
and low limits are the corresponding electron diffusion values or the 
quasi-neutral ambipolar values respectively. 

For the transverse components, we obtain 



D 






- 'X 






-i- 

V -f- pt t- T^a. D 



/ ■*- /.{ ///- s ^ 
+■ D.S 



B 



^—^co > 



(2-li6a) 



(2-!i8b) 



(2-Uc) 






(2-li8d) 
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i 




The maximum transverse internal current will flow vrhen Dm^-p 



are a 



maximum. For the quasi-neutral plasma, equation (2-l;8b), this 
occurs vihen the magnetic field has a magnitude given by 



\ 




For helium, using data from reference [_lT\ , this is approximately 0.12 
veher/n^. The magnitude of the current density resulting from this 
field can be determined from equation (2-22) 



Tr =c(P -r,) ( 2 - 50 ) 

vtiere is a unit vector in the transverse direction. The magnetic 
field generated by this current will oppose the applied field, and 
accounts for the observed diamagnetism of the plasma. This is the 
term ^^^lich satisfies eqxration (2-17), and it must be small compared to 
B to justify neglecting its effects on the other plasma equations. 

In the quasi-neutral plasma, with 7/^ ^71^ , the diffusion tensor 
reduces to 



Oryi-H — 







o 



O 

o Da, 



(2-:;i) 



D. The Continuity Equation 



By substituting equation (2-i;2) into the continuity equation (2-6), 



a second order partial differential equation. For the steady state^ 
this reduces to 



= '^> 1 - ( 2 ^ 53 ) 

whereas for a decaying plasma (for vhich U « 0), it becomes the 

I ' 

diffusion equation 

^ =0 (2-$h.) 

In order to simpliiy equation (2-5h), ve must examine its com“ 
ponents. The components of the diffusion tensor, equations (2-li5), 
are constants because of our assumptions and consequently the 
divergence of the diffusion tensor is zero. In the absence of a 
magnetic field, Dj^+is the scalar D^, and equation (2-52) becomes the 
familiar continuity equation; 



Ds'7'’n. 



Si 



( 2 - 55 ) 



For the steacfy state, ve can further simplify this by employing the 
assumption of proportionality, equation (2-26), to obtain 



V^n . 

17 

In the presence of a magnetic field, ve shall limit ourselves to 



D< 



(2-56) 



orthogonal geometries in which equation (2-53) can be separated into 



DrA,Vl'’K). -t- (2-57) 

In equation (2-57) we have made explicit use of the requirement that 
the transverse density gradient be zero. By employing the assump- 
tion of proportionality (equation (2-26)) and a theorem from vector 
analysis, we see that 



( 2 “!^ 8 ) 



Y\- n " A n / *^ \ K) / ) 

vith a similar equation for the perpendicular components. Thus ve 

can divide equation (2“57) by n_ and expand it to obtain 



If we then define a vector u** such that 





( 2 » 59 ) 




we obtain 

'b* Dmji Ml ^~0 (2»6l) 

which is a first order partial differential equation in u. In the 
absence of a magnetic field, D is a scalar (D ), and equation (2-56) 

D1 S 

can be similarly expressed (in any coordinate system) as 



V^L( 



-h U 



Ji, 

Os 



- O. 



(2-62) 



The significance of the transformation of this equation from a 
second order differential equation in particle density, n, to a 
first order differential equation in u will become apparent in the 
discussion of boundary conditions in a later section . 
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E. Boundary Conditions 
1. Derivation 

We shall now develop an equation for the boundary condition 
and evaluate it for a fevr simplified cases. We consider a quasi- 
neutral stationary plasma in contact with a plane rail. The mag- 
netic field, if any, is parallel to the wall. 



is that the electron and ion currents to the insulator surface 
must be equal; the other is that the current must be continuous at 
the limit of the diffusion region, which can be at the edge of the 
sheath or at the edge of the free-fall zone. 

We have to distinguish different sheath theories accord- 
ing to whether D_, > D+j_ , D.j^ < or D-i ^ , where 



a. The case D^,j_ , ion current inertia-limited. 



Let us first examine the problem in i.rhich D_j_ » 



In this case a negative rail charge builds up until the electron 
current is reduced to equality with the positive ion current. 

The electrons will be treated as having a Maxvfellian distribution 



There are two conditions which we must fulfill. The first 





( 2 - 66 ) 



n . - n^e 




( 2 - 67 ) 
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II 



vhere the potential V at the sheath edge is taken as zero, -e is the 
charge of the electron, and the subscript "o" means evaluation at the 
sheath edge. This notation vill be used throughout this section. 



Ve find that tvro possibilities exist for the ion motion through 
the sheaths (1) their motion may be inert ia-limi ted, or (2) their 
motion may be mobility-limited. In the inertia-limited case, the 
sheath is a free-fall region and the current in the sheath is that 
given by kinetic theorys s 



¥e must next evaluate the diffusion current to the sheath 
edge, which is 



If the electrons and ions are in equilibrium within the body of the 
plasma 



In order that the number of particles varried to the sheath by 
the above diffusion current equal the number carried to the wall as 
expressed hy equation (2-68), the normal component of the diffusion 
current must equal the sheath current. Therefore, we have the 




( 2 - 68 ) 




( 2 - 69 ) 




(2“69a) 



equation 




( 2 - 70 ) 






From this equation we see that 





( 2 - 71 ) 



a result obtained by Ecker [_ 3 ], where, for this case 
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(2-70a) 



oC = ^ 0 ^ ^ ' 

It i-'ill be seen that the boundary condition can always be ex- 
pressed in the form (2-71), vith the value of c<. varying vith the 
sheath characteristics. Let us investigate the potential function V 
for this condition. If there is no magnetic field, ve see from the 
conservation of energy that 

Mv;," _ Mv/, eV 

Z ' -L ^ 

where Mis the ionic mass. From this ve can derive the following 
equation for the ion current density: 



^ KevJ / ^ 

Let us now employ a coordinate system in which the x-direction 
is normal to the sheath and positive towards the wall, and whose 
origin is at the sheath edge. Solving equations (2-67) and (2-73) 
for n_ and n+ respectively, and substituting in Poisson's eq’uation 
(2-l5), we find 



cl V 



e_ 

Zo 



V 






Multiplying by ^ and integrating yields: 






I eV 

kTr 



(2-7J4) 



VThere we have set = -i;,and = 0. This result has been 

reported by Ecker[SJ« For the sheath to be stable, - £ must 

increase with x. Therefore the quantity in the brackets of equation 
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J., 





(2-75) Tiust be positive. Expanding this quantity in series, ve find 






zT^ i x§y 

-r_ 2. tr. 



from vtiich ve see that 




(2-76) 



T; ^ [T_ (2-77) 
is the condition for the sheath to be stable. Bohm {J"} has calculated 
the function V for tvro values of by numerical integration. 

If ve neglect the electron density in Poisson's equation, and 
integrate, ve obtain the Langmuir relation C83 ; 






% 



(2-78) 

b. The case D-j_» D->-^ , ion current mobility- United 

¥e shall next examine the case in vhich the ion current 

through the sheath is mobility limited; i.e., the ions move in the 

field of the negative vail charge according to their mobility. Nov 

the electron current density striking the vail is obtained from kinetic 

theory and the electron density relation (2-67); 

^ eV{d^ 

YioV- e 



L-r_ 



vhere V(d_) is evaluated at the edge of the electron free-flight zone 
(suitably corrected to account for curvature of the electron paths in 
the magnetic field). Again equating this current vith the diffusion 
current brought to the sheath, ve obtain; 

eVa^) 

Vnn 

‘i- / + /-u/u- . 



y)c y- e 



(2-80) 
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In order to further simplify this relation for the case in which the 
electrons and ions are in equilibrium, we use the relation 



from which we see that 




Thus we obtain for the parameter : 



cL - 





( 2 - 81 ) 



( 2 - 82 ) 



(2-83) 



In order to evaluate this expression, it is necessary to know the 
potential function V(x). The ion current density in the sheath 
(neglecting the magnetic field) is; 



y = f X = - /^J.yh e ^ ^ ( 2-8L) 

where 

M ^ (2-81ia) 

Again employing Poisson’s equation, and solving equations (2-8L) and 
(2-6?) for n+ and n- respectively, ve obtain; 



11 . 



e f ^ 
f 4- 




4- 




(2-85a) 

(2-85b) 



At the sheath edge (v4iere x = 0, V = 0, and 

a result consistent vrith the assumption of quasi-neutrality at the 
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sheath boundary. Equation (2-85b) can be integrated once, and yields 
a first order differential equation for V; 




( 2 - 86 ) 



Further integration, to obtain V(x), must be performed n'jmerically. 
Hovever, if we again neglect the electron density, equation (2-86) 
can be integrated to yield; 




(2-87) 



c. The case D+i 

An approximation for the case ^ D+j_ can be easily 
obtained for inertia-limited ion current. By neglecting 
as compared to 1 in the equality 



p- », 

v'e obtain 



( 2 - 88 ) 






from which it is seen that 






B^iploying this relation in equation (2-?0a) yields 



(2-89) 



( 2 - 90 ) 



If vre neglect the last term in the parenthesis as compared vrith one. 
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( 2 - 92 ) 



and employ Einstein's relation 

Jl: - hZi 

' e , 

we obtain 




( 2 - 93 ) 



do The case 

For the plasma in which D+j>> , a positive vrall charge 

must build up until the ion and electron currents are equal. The ions, 
assuming Maxwellian distribution, will have a density vrithin the sheath 
given by 

( 2-910 



where the ionic charge is +e. 

The ion current to the vail is obtained from kinetic theory as 



yicV^e 

~k ^ - 



k n 



‘-f ■> 

v’here the evaluation of V(d+) is at an ionic mean-free path, d,. , 



(2-95) 



vrhich has been adjusted to include the effects of the magnetic field. 
Equating this to the diffusion current to the sheath yields the 
boundary condition (2-71) in vhich 



3 s 

= je )>% 

In order to evaluate the function V, it is necessary to employ 
the equation for the electron current (neglecting the effects of 
the magnetic field) s 
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(2-97) 



By the sane procedure used for the nobility-limi ted ions vhen 
D-j>> D+j_, we obtain 




(2-98) 



2. Evaluation of oC 

The over-all problem of evaluating oC of the boundary condi- 
tion equation (2-71) is a nuch more complicated problem than indica- 
ted ty the derivations of the previous section. The influence of the 
transition region, the effect of non-planar geometries, deviation 
from thermal equilibrium within the plasma proper (and in the transi- 
tion region), and secondary emission are just a few of the factors 
which can affect this quantity. The transition region, trhich joins 
in a continuous fashion the sheath and the m.ain body of the plasma, 
presents particularly difficult problems which are beyond the scope 
of this thesis. In the computation of ^ for a surface i.7hose dimen- 
sions are comparable to or less than the mean free path, it is 
especially important to compute the fraction of particles brought 
to the sheath edge which would miss the surface in their flight 
through the sheath and would pass back into the plasma. Thus, for 
the limit where the surface is very nuch smaller than the mean free 
path, as for a Langmuir probe, can become vanishingly small. 



3. Discussion 

An understanding of the processes occuring within the sheath 
has been shovm to be necessary for the evaluation of the boundary 
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condition equation (2-71). Ve have seen that the requirement that 
the current be continuous results in this condition i-hich must be 
met at the surface of the ion free-fall zone (the region vrithin one 
ionic mean free path, , of the vail) or at the surface of the 
sheath, whichever is farther from the vail. As the sheath thickness 
and X*- become very much smaller than other dimensions, equation 
(2-71) nay be replaced ty its limit for vanishing /S+ : 



This is the generally-used boundary condition, but for small values 
of and large , the error in (2-99) can be appreciable. From 
(2-71) we obtain? 



The importance of equation (2-71) and of the transformation of 
the continuity equations to the forms (2-6l) and (2-62) is nov 
readily seen. The normal component of the vector u defined ty 



edge. If there is more than one vail, and the vails involved are 
orthogonal, the density gradients established to satisfy this condition 
at each sheath edge vrill not interfere with each other. Parallel 
surfaces, on the other hand, create an eigenvalue problem because u 
must satisfy the boundary condition at one sheath edge, must change 




(2-99) 




( 2 - 100 ) 



where the subscript refers to evaluation at a distance A^from the 




2h 



as prescribed by the continuity equation (2“6l) as ve move across 
the plasma, and must acain satisfy uniquely the boundary condition 
at the other surface« If the walls are neither orthogonal nor parallel, 
the density gradients established to meet the boundary condition for 
one surface will have along the other surface normal components which 
are not constant, and vice versa. This will result in a violation 
of equation (2=71) » Consequently, for non-orthogonal geom.etrics we 
are unable to predict the steady state distributions with the equations 
developed in this thesis except ty an approximate method to be 
discussed in a later section, 

F. Heat Losses 

The only source of heat loss considered in this thesis is that which 
occurs at an insulated surface as a consequence of the particles 
impinging on the surface and giving up their energy. This energy can 
come from several sources? (1) the recombination energy, eV£, from 
each electron-ion recombination, (2) the kinetic energy of the particles, 
and (3) excitation energy of excited ions which return to their stable 
state upon colliding with a vail. All of these are proportional 
to the particle current density to the surfaces. As a simplification 
for illustrative purposes in later calculations we shall limit ourselves 
to plasmas whose temperatures are sufficiently lov' (belov* about 10^ 
degrees) that only the recombination and electron kinetic energies 
are appreciable. The heat liberated at a surface is then 




( 2 - 101 ) 



v'here the subscript "n” means normal to the surface and is 

expressed in electron volts. 
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SOLUTIONS OF THE CONTINUITY EQUATION 
A. General Solutions 

This section discusses a procedure to obtain particle distri^ 
bution functions in several coordinate stems, and presents some 
of the general features of the solutions. The follovang sections 
develop many solutions for specific geometries. 

1. Exact Solutions 

a. General Development s 

With no magnetic field, vre have exact solutions in the 
eleven orthogonal coordinate systems in vhich the vrave equation is 
separable vthen the boundary surfaces coincide vrith coordinate planes. 
The procedure ve shall describe applies to the non-steady state 
plasma for vhich we repeat the general continuity equation (2-^5) 



which becomes a wave equation for the steac^-state 



±L 

for a decaying plasma vrith no ionization it is 



D.V%. = 1 = , 

with a steady-state counterpart (Laplace's equation) 



(3-1) 



(3-ia) 



(3-lb) 



V 



2 

71 - 



= O. 



(3-lc) 



In the presence of a magnetic field, ve further require that the 
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coordinate system be one in vhich the continuity equation may be ex« 
panded as 



>1. -b 4- (3-2) 

St 

(see equation (2“57), 

Eqtuations (3=1) and (3“2) are separable in the geometrical coord- 
inates and time, and the general solution for a generalized coordinate 
systems u, v, w, may be expressed as 

n(u,v,w,t) s N(u,v,w,)T(t)^ (3-3) 

S’obstitution of equation (3=3) in equation (3“2), and division by NT 
yields 



^ A/ ^ i: , 

N -~FT' * ^ ' r = 



(3~ii) 



vhere is a separation constant. The time solution is 







(3-5) 



and then 



D«„V„"W 

" —jr- 

In the absence of a magnetic field, equation (3-6) reduces to 




in idiich 






(3-6) 



\e fiJ 
lz 

X = Ds 

^ / 

and, for the steady state (kt ® O) 






1?^ - 1. 



S t> 



(3-7a) 

(3-7b) 
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Our approach to the solution of equations (3“6) and (3"7) is to 
employ the technique of separation of variables. Thus we assume 
that the solution is of the form 

N(us,v,v) ® f j(u)f2(v)f^(w)j (3=8) 

and, by substitution into equation (3”6) or (3“7) (as appropriate), 
we obtain three differential equations, one in terms of each of the 
coordinate variables, which the functions fj, f 2 > and f 3 must separately 
satisfy. Magnus and Oberhettinger 1^93 have summarized the results for 
eight coordinate ^sterns. Because the separated equations are of 
the second order, there will be two linearly independent sets of 
eigenfunctions for each coordinate, F^^^(k^,u) and ,u), 

G^^I%i,v) and G(2)(kmjv), and kn^w) and H^^^(kn,v), which satisfy 

the differential equation in the corresponding variable. We shall call 
the functions with the superscript (1) functions of the first kind 
and the othei s functions of the second kind. The most general solu« 
tion to these equations will then be a sum of the products of all 
such linearly independent solutions of the form 

X [Ck, H" (K H 

The complete solution for the particle density distribution is 
then 



(^- 10 ) 

(3-iOa) 
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The summation is over all possible values of m, and n, k^, 
and kf, are eigenvalues or spearation constants (depending on the 
coordinate system chosen) corresponding to m^ and n respectivelyo 
¥e can sum over only those eigenfunctions vrhich belong to a degenerated 
eigenvalue (as determined try the boundary conditions). It should be 
noted that the time=dependent solutions to be developed herein are 
for those problems for the boundary condition (2=71) or (2~99) 

is still applicable. A plasma decaying from a steady-state distri- 
bution upon cessation of ionization is considered to be such a 
plasma. 

b. Cartesian Coordinates 

An example of the above method of solution is the solution 
in cartesian coordinates in thich the functions of the first and 
second kinds may all be cosines and sines respectively j or in vhich 
one or tvo pairs of the functions may be exponentials. For example? 

^ (3»n) 



By substitution of this solution in equation (3=7) j vre see that 
k^ of equations (3~7a) and (3"7b) is given by 



and consequently that the time solution (equation (3=5)) is 

-r , ^ 
iP'ntv, “ he. ■ 



(5=lla) 



(3=llb) 
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Obvious lys when 




( 3 - 12 ) 



the tine dependence vanishes^ and the plasna is "stationary«" 
Another possible solution^ involving exponentials, is 




where is given by ^ 



(3=13a) 



c. Cylindrical Coordinates 

There is an independent eigenvalue, kj, for each set of 
functions only in cartesian coordinates. For any other coordinate 
system u,v,w, defined by a transformation equation of the form 



those eigenfunctions of the variables u,v,v which appear together in 
ar^ of equations (3-llt) will be connected through a separation variable 
which is required to reduce equation (3“ii) or (3-7) to separate 
equations in each of the coordinates. Each such separation variable 
required for the above purpose reduces the number of eigenvalues 
available for the satisfaction of the boundary conditions. For cx= 
ample, for cylindrical coordinates, where 



the eigenfunctions in r and are related by the separation constant, 
m, as indicated in the xfcll=known solution 



X ® x(u,v,w) 
y ® y(u,v,w) 

Z z(u,V,V?)j 



(3-lli) 



X ® rcooif 

y ® r 



z = z 



(3-15) 



30 



( 3“16) 



A/Co 'C ar) k) 



Trtiere 









and 



= C^^dc3a(7vifKC5Vi^6^^') 
- Cy 4? C<- ^ /, V 

2 ('«;) = Ciy^CcaClc'n^') ^ Om 



^ ' ^v‘ ^ ^ ( 3-l6a) 



M i 0 

yV) • O 

l^n'^ ^) 

fcy,^ O 



'y 



- f 6/t ci 

An alternative Z(z) solution is 



( 3“l6b) 



t( e) ‘ ^ ^ ^7yt 



K'^0 



(3- 16c) 



vith 



K\ (3-161) 

J,^( k r) and l^(k r) are Bessel functions of order m of the first and 
second kinds respectively (Bessel and Neumann functions), 

c. Spherical, Prolate and Oblate Ellipsoidal Coordinates 

In the other coordinate systems to be discussed, the general 
solutions of (3“7) ares 
Spherical s 



R(r)Q(e)i(f\ 



(3-17) 
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vtiere 



R(f) ‘i-'" 



= G,r'’ 4- 



Ifj - 0 



(S)(b) = pycu:>0) 

^ ' Q -t (/) ho /a 

= Cc,^^Cj 



*vi r O, 



(3-175.) 



and 



k^-.pA 

and in tdiich P^(cos©) and Qn(cos©) are Legendre Polynomials of the first 
and second kinds respectively and jjj(k r) and nj^(k r) are spherical 
Bessel functions of order n of the first and second kinds respectively® 
The latter are defined fcy the equations 

Prolate Ellipsoidal (See Appendix I for diagram and transformation 
equationsi notation employed is that of Flammer m vith slight modi- 
fication) s 

(f, t Cn„„ %nti 

vhiere 



(3-!?c) 



and d is the interfocal distance, and 4n^(c^, are sphero i da I 



32 



angle functions of the first and second kinds and ^ 

Cj^ ^ ) are spheroidal radial functions. See Appendix II for 

tables and graphs of these functions. 

Oblate Ellipsoidal functions nay be obtained fro-n the prolate func- 
tions by the substitution of =ic^ for where appearing in equation 
C3-l8). Equation (3-l8a) is unchanged by the transformation. 

Similar equations are obtained in parabolic cylinder coordinates 
(where the functions include Weber-Hemite functions), parabolic 
coordinates (V-hittaker furctions), and elliptic cylinder coordinates 
(Mathieu functions). These are not developed further in this thesis 
because their geometries are not useful in the problems with vhich ve 
shall deal. 

d. Solutions in the Presence of a Magnetic Field 

In the presence of a magnetic field. Dm is a tensor, and 
equation (3-6) must be employed. In order to make a separation of 
the form (3-6), we are limited to cartesian coordinates or cylindrical 
coordinates for the problems lAiich we shall treat. Parabolic cylinder 
and elliptic cylinder coord iiiates are also suitable, but we shall not 
need them for the geometries with which we shall be concerned, Fv-rther, 
if we consider the assumptions under ^ich Dn vras developed, we may 
immediately limit the treatment to the case ^Jhere m and C are zero, 
since no gradients in the ^ direction can exist. 

In the coordinate systems selected, equation (3-6) can 
be reduced by separation of variables, by assuming 




(3-19) 
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and thus obtaining two equations 



and 








where 



( 3-20a) 



(3-20b) 



The solution to the above equations is the sane as in the absence of a 
field, except that equations (3-12) and (3-l6d) become (for a field in 
the z~direction) s 



k \ (3^12«) 

and 

^ f I ? I (3-l6d«) 

e. Further Simplification of Previous Solutions 

The equations developed so far can be simplified considerably 
considering some of the restrictions upon the solutions* The require- 
ment of a single-valuedness limits/, n, and n to be integers, and it 
is possible to further limit then to positive values without loss of 
generality. Further, if the plasma possesses rotational symmetry 
(vdiich it must in the presence of a magnetic field under our assump- 
tions), then m and C5 are both zero in the equations for N(u,v,w) 



(vith the exception of rectangular coordinates). The requirement 
that the density be non-negative everywhere within a given region will 
further limit the solutions. Because all the functions discussed 
(except the exponential function) are periodic, this latter requirement 
usually limits the stea<^-state solutions for bounded regions to the 
first positive segment of the applicable function, and values of , 
kj,j, and/or are thus determined by the dimensions of the plasma 
container. From the boundary condition (2-71)j the physical reason for 
this effect is apparent. u =— has a definite value at the edge of 
the sheath. Equation (3“?) can be transfor^ied into an equation in u 
in the same manner as in the derivation of equation (2-62), yielding 

V-U (3-21) 

It is the value of k^ in this equation that determines the "trajec- 
tory” of the vector u across the plasma to the sheath cn the "opposite" 
side where u must again satisfy the boundary condition (2-71). Because 
equation (3-22) is a first order equation inli, the problem is completely 
stated (except for a condition to determine an arbitrary normalization 
constant), and only one value of k is possible. 

Further, and obvious, simplification can be accomplished for many 
of the solutions stated herein if the>'' must describe a plasma in a 
region in ihich some of the eigenfunctions are sincular. In parti- 
cular, the following functions have singularities for the values indicated 

QS(ii), s<9(o,ii), s<n^ (-icUi). 
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2. Approximate Solutions 

There are a limited number of geometries for vhich a stable 
ambipolar solution may be obtained by the separation procedure de- 
scribed herein. It is, however, possible to obtain an approximate 
solution for the particle density distribution if the geometry is 
not too radically different from a geometry for which an exact 
solution is ciitainableo 

In the discussion of the boundary conditions (Section II-E-3), 
it was shown that the conditions at the sheath edge control the den- 
sity distribution in the vicinity of a surface. If there are tvo or 
more surfaces bounding the plasma in such a manner that it is not 
possible to obtain an exact solution for the particle density dis- 
tribution, then ve divide the plasma into "zones of influence" 
about each of these surfaces. ¥e then calculate the density dis- 
trihition for each zone as determined ty the boundary conditions at 
the insulator surfaces belonging to that zone. In the region near 
such surfaces we expect our solution to be very nearly correct. There 
will be a transition region between the zones of influence in vhich 
we expect the greatest difference between the model and the actual 
plasma. In this region there is either a departure from the ambi- 
polar diffusion or a time dependence (non-stat ionary plasma), or 
both. However, we expect that processes in this transition 
region do not affect the total solutions significantly in many 
problems. A measure of the importance 
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of these processes is the degree to vhich the solutions for the 
different zones conflict in the transition region. When the dis- 
continuity across the transition region is small, we expect that our 
solution correctly describes the plasma in the immediate vicinity 
of the insulator walls and describes it approximately (vrithin an 
order of magnitude) even in the transition region. 

In order to discuss this approximation more quantitatively, let 
us consider two walls, Wj^ and ¥ 2 , bounding a partly- ionized gas, and 
assume that the configurations are such that the model described in 
this thesis can be applied to each wall in the absence of the other. 
Thus, two distribution functions, nj(qj) and n 2 (qj), arc obtained 
(with qj being the coordinates). If, through adjustment of the 
spearation variables of equations (3-6) and (3-7) and of the normali- 
zation constants, it is possible to make 

(3-2lia) 

and 

everywhere, then the solution is the exact solution. If this cannot 
be done, then we must find a surface, S, which separates the zones 
of influence of the walls and ¥ 2 . The separate solutions for 
these zones, joined at S, form the zero-order approximation of our 
solution. 

There are three general procedures ty vhich we might develop a 
criterion for determining the position of the surface S. ¥e shall 
briefly discuss all three in the following paragraphs, but shall not 
attempt a rigorous proof. 
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One procedure is to employ a perturbation theory in which we 



define 



and 


77^ ^ rio-t ^ ' ' 


(3-25) 




£ = Eo - 


(3-26) 


where ){ 


is a parameter defined hy the equation 






77^,V'i -- 


(3-27) 




in which 


n^Q is defined as the maximum density in the plasma. 


From 


equation 


(3-27) and Poisson's equation it is apparent that 

\J ^ o 


(3-28) 




eXo 


For 


our zero approximation, we impose the conditions; 




and 


7-P^-pk = o 


(3-29) 




r„ 


(3-30) 



in each zone of influence. 

In order to determine a criterion for positioning the surface 
between the zones of influence, S, (i,e,, in order to find an adapted 
Zero-order approximation), we must look at the first order terms. 
Neglecting quadratic terms in , we obtain 

1 

f 

+ 7?o F 'l - D+Vt?+ (3"31a) 

and 

P,' . -/U.(£oy>- D-Vt).' (3-31b) 

[ 

i 
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If we further assume that either there is no time dependence or that 
the time dependence can be accounted for by replacing all quantities 
in the following equations ty their time-averages, we obtain the 
following from the continuity equation: 

V'Tl' -hm't (3-32a) 

and 

V-f.Vp^y . -aH) P-4') (3-32b) 

where the subscript ”S" means evaluation at the surface S, 
is the Dirac delta function, and n is the unit normal to the surface 
S. The last two terms in the above equations represent particles 
which wuld have to be created or absorbed at the imaginary surface 
S in order to supply the currents indicated. The first order density 
terms can be either positive or negative. It is apparent from (3-32) 
that at the surface S, 

jlPo,, ^ -44)= -(P:,. -K) (3-33a) 

If the terms of equation (3-33) are identically zero, then the 
first order solution can be the trivial solution 

r?t =0 (3-3lr) 

Any deviation of the first order solution from equation (3~3h) is 
obviously a function of the magnitude of the first term in equation 
(3-33). The surface S is thus the surface at v'hich this term is a 
minimum. If there were a surface along which 

<3-35) 
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everywhere, subject to the condition that the density solutions were 
also equal at the surface, then the ^ would 

vanish. This is not possible generally, however, and a least-mean- 
square procedure should be applied. Because this is not feasible 
for many problems, the following criterion for determining the 
position of S, applicable at a given point in which we are interested, 
is employed to obtain our adjusted zero-order approximation; 



iDrw 




(3-36) 



where / 1 is the angle between Vf\o\ and the surface normal and iT'i 
between and the surface normal. Further, if 







(3-37a) 



and 



then at the surface S 

Voil ^rioz^ -C-f, ^ 

Solving for C 2 , substituting the result in equation (3“36), and 
dividing both sides by Cifi(qi)s yisi^s 



(3-37b) 



(3-38) 











(3-39) 



Because only order of magnitude results will be obtainable from 
this approximation, it is reasonable to assume as a first approxima- 
tion that 






(3-iiO) 



ho 









3 ? 



By employing equations (3"37)> V'e obtain the follovdng conditions 
defining the surface of separation 

I 'p. 

I Uno -yi^^ 

The second procedure to derive equation (3-Iil) is to consider 
the relative influence of each boundary at all points v’ithin the 
plasma, and define the surface S as that surface at vhich the in- 
fluences of the two boundaries are equal. Some measure of this 
influence must first be found. If we have in the neighborhood of 
a boundary a solution for the particle density vrhich has a maximum, 
then certainly beyond that maximum the plasma is not influenced 
appreciably by that boundary. For the region between the surface 
and the maximum, the measure of a boundary's influence could be 
the current created in response to its presence. This procedure 
leads to equation (3"lrl) with an added restriction that if there 
are two points at vhich this equation is satisfied, v^e must select 
the one within the zones defined ty the walls and the first maxima 
of their respective solutions vhenever possible. 

The third procedure to determine the position of the surface 
S would be to calculate the probability that a particle at a certain 
position would go to the one wall or the other. The surface S 
dividing the two zones would be the surface along which particles 
would have equal probabilities of reaching all the walls in question. 
This might be calculated by stochastic methods or from the diffusion 
equations. ¥e did not attempt the former. The latter method is 
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(3-hl) 
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to find the surface from which the mathematical solution indicates 



that as many of the particles produced vrithin a 
straddling the surface would go to the one zone 
This procedure also leads to equation (3-i4l) as 



differential voltime 
as to the other, 
the criterion. 
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B. Specific Solutions of the Continuity Equations * 

1, General Conditions 

In order to obtain an estimate of the magnitude of certain 
quantities in the equations to follow, we shall use the follovdng 
sample data, some of which is taken from 3roi.m J 
gas; Helium 
p s 1 mm Hg. 

T « 700° K 

T = 50,000°K eV, 



ps 10^ cm^/volt-sec js 1 m^/volt-sec 

7 o 6 X 10 ^ cm^/volt-sec — 76 m^/volt-sec 
D^. m^/sec 

D_ 5 s ^ ^ 3U0 n^/sec 

Ds ^ Dg = P-th. = U.li n^/sec 



D«- D+ 

eVi ^ 
eVj + 



^ Uoli volts 
2lio^8 eV 
* 29 eV 



E ® 300 volts/m 



('3-50) 









” Hq “ 10^^ cni~^ 



1o18„-3 



where '»r^’ is the density in the center of the plasma or in an un“ 
disturbed, constant-=density region. 

For the calculations of this section \jz assume that the mean free 
paths of the ions and electrons are both sufficiently smaller than ary 
boundary dimensions involved, so that boundary condition (2-99), 






is applicable. The more general boundary condition (2-71) will be 
employed in the discussion of cylindrical probes in the next chapter^ 
and the method of employment of that boundary condition will be dis- 
cussed then. 

A quantity often quoted and employed in the literature is the 
diffusion length, A, i-hich is defined by the equation 




v.z shall compute this quantity for the various geometries in this 
section. 

Finally, ’■^e shall calculate the "zone of influence" of many of 
the surfaces to be discussed. This we define as that region in i-hich 
the particle density is less than 50^ or less than 90 % (as indicated) 
of the density at the center or undisturbed region of the plasma. 

2. Plane and Rectangular Configurations 
a. Rectangular cavity 

Because it is the most complex geometry we can treat in 
rectangular coordinates, we begin vdth a discussion of a rectangular 
cavity of dimensions Xq, yo, and Zq. The limitations of the develop- 
ment of the diffusion tensor prevent the application of the follov?- 
ing solutions to magnetic fields, and thus the scalar coefficient Dg 
vrill be used. Of the possible variations of the solution (such as 
equations (3-11) and (3-13)), we can immediately rule out those con- 
taining the exponentials because of the requirement of zero density 
at the walls. The boundary conditions also limit the solution to 
the first positive segments of the sine functions (and exclude the 






cosine functions) if the origin is taken at one corner of the cavity. 

If the origin were selected at the center of the cavity, the cosine 
functions would be used and the sine functions excluded. V/ith the 
origin at one corner, and the sides along the positive coordinate axes, 
we obtain the follovring solution for the steady state = O) plasma; 

^ (3-52) 



with the diffusion length given by 




n 



^ JL ^ (TLY 

Ds ' 
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^ *(fy 



(3-52a) 



The three components of the particle current density at any point 
are easily obtained, using equations (2~h2) 

Ds 7; 

with similar equations for the other components. From these equations 
it is seen that the particle current density is zero at the center of 
the cavity, is directed outward from the center to the vralls, and is 
a maximum at the vralls where, for the data of (3-50) and a 5 centi- 
meter wall-to-vrall separation, the magnitude of the current density 
at the center of a wall would be approximately 2.7 X 10^^ particles/cra^ 
-sec. It is of interest to note that the steady -state particle current 
densities at the walls decrease in inverse proportion to a linear in- 
crease in cavity dimensions. This decrease is a consequence of the 
decrease in the gradient of the current density. From equation (3-52a) 
one can see that the rate of ionization, )) , necessary to replace the 

particle losses to the walls is inversely proportional to the square 
of the dimensions (and, incidentally, is primarily determined by the 



smallest dimension of 



the cavity as a consequence) » The explanation for the inverse 
dependence of the rate of ionization per unit volume is easily seen 
by equating the number of particles destroyed at the vails to the 
total number created vithin the cavity. The latter is proportional 
to the volume and to y , The number of particles carried to the vallf 
hovever, is proportional to the surface of the vails and to the current 
density. The latter is inversely proportional to the linear dimen- 
sion, as previously mentioned, and consequently the number of parti- 
cles drawn off is proportional only to the first povrer of the linear 
dimension of the cavity. Consequently, the inverse dependent of 
on the linear dimensions is expected. The fact that the number of 
particles destroyed at the vails is directly proportional to the 
linear dimensions of the cavity (and therefore, so are the energy 
losses) is also worth noting. 

The space charge field, according to equation (2-36) is 




or, again using the data from (3"50), is approximately -270 volts/m 
(directed away from the center of the cavity) at the point (.2 ^Xq, 
,5yo> o^Zq). The cavity dimensions are assumed to be 5cm. Again, 
an inverse relationship to the wall-to-iv-all separation exists. The 
space charge potential (equation (2-1^1)) is 

or approxiinately =1,53 volts at the point (.25x^, .^y^, ®5 Zq) for 
a5 cm. cavity. The potential is with respect to the center of the 
cavity. The potential predicted hy equation (3'=5ii) becomes singular 

U6 



at the vails, and this feature of the solution indicates a limitation 
of these equations? they are not applicable beyond the sheath edge. 

Using the current equations (3"53) we can compute the energy 
lost at the walls. From equation (2-101) ve obtain 












(3-55) 



viiere the factor 2 takes into account the surfaces adjoining the corner 
opposite the origin. This equation reduces to 









(3-55a) 



7Y ^ Xu'jo^«> 

The heat liberated per unit area is obtained ly dividing each integral 

in (3-55) by the appropriate area, obtaining for the three sides of 

the cavity: 
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(3-55b) 



Again using 5 cm, dimensions and the data from (3-50), vre see that 
^’all approximately i*,86 X 10^^ e'^/sec js- 7»B joules/sec. From 
equation (3-55a) it is seen again that the total heat lost per seend 
is proportional to the linear dimensions, and that doubling the 
length of all the sides, for instance, x^ill double the rate of heat 
loss. 



One very important fact about plasma characteristics nay be seen 
from the preceeding equations, and in particular from equation (3-52a). 
That is tliat the effect of a small dimension on the characteristics 
is very much greater than that of large dimensions, so that a long, 
narrow plasma has, in the middle region of the longer dim.ension. 



characteristics that are almost completely dependent on the shorter 
dimension and are relatively independent of the longer dimension, ThuSj 
a plasma t-hose z-dimension is very much greater than its dimensions 
in other directions can often be treated quite properly as ''infinite'* 
in the z=directicn , with corresponding simplifications of the equa“ 
tions describing the plasma. 

The zones of influence of the rectangular surfaces are indicated 
in Figure 3=1. 

If ionization ceases, P ■e 0, and the plasma decays exponentially 
in accordance vith equation (3=^), with given by (3“7a), (3-12), 
and (3“52a) as 



t.’hich, for a cavity of 5 cm dimensions, is about 20 ^sec. This 
decay time is inversely proportional to the square of the linear 
dimensions. 

If the iritial distribution is not the steady-state distribi- 
tior, but instead some function f(x,y,z), then the coefficients of 
equation (3=11) must be obtained in the normal manner for a Fourier 
series, and a series of the form (3=10) must be developed. For example, 
if the initial distribution is a constant density throughout 
(except discontinuously zero at the surfaces), the particle distribution 
function will be 




( 3 - 56 ) 



Thus the plasma will decay to l/e of its initial value in time. 




( 3 - 57 ) 







( 3 - 58 ) 
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with 




(3-^8a) 



It is to be noted that the higher order terms decay at a faster rate, 
so that the plasma approaches the steady-state distribution as it decays. 
Figure 3-2 shows the decay of such a plasma at three different times. 
Other features of this solution are that the rate of decay is in- 
versely proportional to the square of the linear dimensions so that 
a plasma in a large cavity decays more slowly than one in a small one. 
Further, vre again see that the shortest dimensions control the decay 
rate, so that dimensions vihich are very much larger than the smallest 

dimension may again be treated as infinite. For instance, for a 
0 

length which is ten times that of the, others, the error in ignoring 
it in equation (3=57) is less than one percent. 



b. Infinite tube, rectangular cross-section 

^y- allowing Zq to go to infinity and setting s ^y., 

in the equations of the previous section, we obtain the equations for 
an infinitely long rectangular geometry. Thus w have for the steady/ 
state 



vith 



(3-59) 




A’- 





(3-59a) 



The current density, potential and other equations are similarly ob- 
tained and vdll not be stated here. 
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c. Infinite Plane 



By again allowing one dimension to become infinite, for 
example, Yq, and by replacing y by Yq/ 2, we obtain the equations for 
the simplest geometry, two parallel planes separated by a distance x^, 
or a single plane in contact with an infinite plasma in which 



where we know the rate of ionization from other considerations, or 
where Xq/ 2 is the distanceto the undisturbed region of the plasma. 

The restriction that there be no magnetic field can now be dropped if 
the only density gradient is along the coordinate lines perpendicular 
to the planes and if the field is parallel to these planes. If the 
field is perpendicular tc the plane (so that it does not affect diffusion 
perpendicular to the plane), or if there is no field, then Dmj^ in 
equation (3-60) and in the equations to follow can be replaced by its 
no-field limit, Ds, 

The particle density distribution for the steac^ state is 






(3-61) 



with 








(3-6la) 



The other equations are derived from those describing the rectangular 
cavity by replacing y and z by Yq/2 and Zq/ 2 respectively, by lettJng 
yo and Zq become infinite, and by replacing Ds by Dmj_ . 



3. Cylindrical Conf i gurat ions 
a. Finite Cylinder , interior 

For a cylinder of radius R and length z^, we can immediately 
simplify Equation (3=’l6) for the steac^ state by considering the 
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functions and the boundary conditions. The Neumann functions are 



singular at the origin, and so their coefficients, , must be 



zero. Further, if the origin is taken at the center of the base of 
the cylinder, the coefficients of the cosine functions > ‘"Tn* must be 
zero also, because the density is zero at z = 0 (using boundary 
condition (2-99)), Then, as in the case of the rectangular configura- 
tions, the requirements that n alvtays be positive (or zero) and the 
boundary condition at z « zq requires us to set all the sine 
coefficients to zero except and becomes . Because we 

have discussed no forces vhich would tend to esteiblish an assymetric 
angular distribution, and indeed in the case of the derivation of 
the diffusion tensor ir a magnetic field have specified that no 
such as^ymetiy exist, we set m = 0. Thus vre are limited to the 
zero order Bessel function, and because of the boundary conditions 
must select the first positive segment of that function. Thus we ob- 
tain the well-knor^m distribution function 




( 3 - 62 ) 



where 



-L 




and 




(3-62b) 



The magnetic field, if any, is parallel to the cylinder axis. 



The ion current density components are obtained ty using 
Equation (2-U2); 
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(3-63) 



Thus, at the center of the cylintier (r = 0, z = ^ z^), the current 

(density is zero, and within the remainder of the cylinder ) and i 

are both directed from the center to ,the walls, | ^ is a maximum 

ri '■ 

at the ends of the cylinder, but 1^.. reaches a maximum (see Figure 
3-it) when .. 









. { 



i,e,, at r s 0.75R* This maximum occurs at a radius less than R 
because the additional current supplied ty ionization b^ond r « 0.75R 
is not sufficient to counteract the effect of the increasing cylin- 
drical volume into vrhich the current is flowing. As for the rectangu- 
lar cavity, an increase in cylinder size is accompanied by a decrease 
in ionization rate and,. current densities in the steady state. 

The space charge fj eld, as obtained from (2-35) and (2=36) is: 

( I ~h/A ~ D - ,2. V 

^ (3-6U) 



Hr ^ £j. 




The zones of influence of the cylinder surface are indicated 
in Figure 3-3* 

One of the characteristics of this plasma vrhich is of interest 
is the fact that the shape of the steacfy state configuration of ion 
density is dependent only upon the cylinder dimensions if the strength 
of the magnetic field is less than that at vhich the "pinch" effect 

5U 



begins to appear. However^ the ion current density in the radial 
direction (perpendicular to the field) is proportional to the inverse 
square of the magnetic field (at high B) through its dependence on 



If the steady state plasma begins to decay, the distribution 
is given by the equation? 

_ (3-65) 



If, however, a decaying plasma vrith initial distribution F(r, ^|^,z) 
is to be described, it is necessary to develop a series from Equation 
(3-l6) in the same manner as used to develop a Fourier series, obtain- 
ing equations similar to (3-^8) , 

b. Infinite Cylinder , interior 

In a manner similar to that of the rectangular problem 
it is possible to treat a very long cylinder as if it were effectively 
of infinite length. Employing identical procedures, we obtain for 
the particle density distribution function of a "long" cylinder of 
radius R? 



7l(r) - 7?o 



where 



V - (^) b«j_. 



-A’ 

The particle current density is given by 



( 3 » 66 ) 

(3“66a) 



P, - + .. J, r) , r, = P, -- O, 



(3-67) 



the field by 
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( 3 - 68 ) 



the space charge potential ty 



V- 



M- 






(3-69) 



and the heat liberated at the vralls per unit cylinder length try 



Eo^cJd/ 

/u/ 



(3-70) 



^ 1 ^ 6^71 

It should be noted that the latter is independent of the cylinder 
dimensions, and using the data of (3-50), is approximately 1.32 XIO^^ 
B"2 eV/m-sec = 2.11 joule/m-sec (with B in vebers/m^). 

The decay problem is treated in the same manner as for the finite 
cylinder, and we obtain for decay from the steacty state 



Figure 3-h is a plot of the radial current density as a function 
of radial position for a steady plasma. Figure 3-5 is a graph of 
the space charge potential versus radial position for an infinite 
cylinder. The fact that the space charge field and potential become 



infinite at the wall is a consequence of the assumption that n^r^ii 
= 0. As in the rectangular case, these equations are not valid within 
the sheath. 



c. Infinite Cylinder in Infinite Plasma 

Using our tectiniques, there is no exact steady-state 
solution to the continuity equation (3-1) or (3-2) for a plasma exterior 
to a cylinder of finite length because we cannot combine the solutions 
(3“l6) in a manner such as to meet the boundary conditions at the ends 
of the cylinder. However, for the infinite cylinder of radius R 
it is possible to solve for the distribution function. Again ve limit 
ourselves to those problems in vAiich no assymetric forces are acting 



f 





0 ^ 





TTITT 







(and thus to axially syitmetric plasmas), but retain the Nevanann 
functions because the origin is no longer within the limits of r. VJe 



obtain from (3-16), for the steacfy state, 






(3-72) 

vrhere A is a nomalization constant which we shall determine later, 

R is the cylinder radius, R“ is that radius at vhich n(r) is a maximum, 

and . 

"V ( i-kL o \ 

(3-72a) 



C = 



if the boundary condition (2-997 applies, and is a more 

complicated function if equation (2-71) applies. The latter function 
will be developed in the section on probe theory, the radius at 
vihich equation (3-72) is a maximum, is found by setting tlie first 
derivative to zero| i.e,, it is the solution of 



j; (7^, R"") - CW, 



(3-72b) 



The normalization constant. A, can be determined if the density n(r|^) 
at a radius rj is known. Then 



A 



'7 ( r, 






The radial current density is 






Ji 

0. 



I 



r 






(3-72c) 



(3-73) 



'I' ■ 

d. Coaxial Cylinders 

The equations of the previous section also describe the 

particle and current density relations which exist between the surfaces 

of coaxial cylinders of infinite length with inner radius R and outer 

radius R' . Because of the boundary conditions, y is no longer an 

independent variable, but is an eigenvalue in the steady state plasma, 
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Using the boundary condition of zero density at the vails, the value 
of is seen from equation (3‘=72) to be that value for vdiich 

+ CWo(vSiR’') = (3-™ 

ii. Spherical Geonetries 

Distribution functions for decaying plasmas with arbitrary 
initial distributions and for steady-state plasmas in spherical geo- 
metries can be obtained hy series development of equations (3-17). 
Equations for the spherically symmetric plasma in the steady state 
(^^»ith ionization occuring) are very similar to (3-66) and (3-72) vith 
spherical Bessel functions (equation (3-17c)) replacing the cylindri- 
cal Bessel functions in equations describing spherical geometries corres- 
ponding to the cylindrical geonetries treated. In addition, it 
is of interest to describe the steady state solution for a steady or 
very slovly decaying plasma in x^hich there is no ionization occuring, 
and in xhich a sphere has been placed. Ye treat the latter problem 
first. 

a. Sphere in Homogeneous Plasma-no ionization occuring 
If ve have a spherically symmetric infinite plasma 
vith an undisturbed particle density, n^, in v^hich //= ^ = o , 
vith neither magnetic nor electric fields imposed, ve have from 
equations (3-17) (vith iy - m - 0) the folloving density distribution 
around a sphere of radius K 
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V'e have employed the boundary conditions that n(r) be zero at the 
sphere and equal n^j at large distances. The zones of influence 
for such a sphere are 

50fo - 

90fo - r^-loR, (3-76) 

and the radial current density is 






r - 



(3-77) 



This current density is a maximum (of magnitude 



■^qDs 

R 



) at the surface 



of the sphere and approaches zero as l/r^ at large r. The current is 
directed torard the sphere and represents a flov of particles from a 
source at infinity to a sink at the sphere. 

The heat liberated at the surface of the sphere is 







DAeVi 
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(3-78) 



It is proportional to the sejuare of the radius of the sphere. 

Using the equations thus stated, it should be possible to calculate 
the particle density distribution which a small spherical surface in 
a large, slowly-decaying plasma vrould establish about itself, 
b. Spherical Shelly-internal plasma 

From equations (3-17) it is also possible to obtain 
easily the steady state relations for a plasma bounded fcy a sphere of 
radius R centered within a spherical shell of radius R* . V^'e shall 
again assume spherical symmetry, thus limiting ourselves to the 
spherical Bessel functions of zero order, 3y employing the follow 
ing relations for these functions 

; a A. 



ylo ( 1 



Wi 
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(3-79) 



and doing some manipulating, we obtain 



. . ^ r TT('r-R'^l 

/!(r) = C, r ^ J 






(3-80) 



where 



ft " --(i*): 

and for vihich derivation the condition of vanishing density at the sur- 
faces T-ras employed, is an arbitrary normalization constant which 
can be determined if the density is knovm at some radius r2^ within 
the plasma. 

Then, 

= (3-60b) 

and 

( 3-800 

The radius of maximum density^ R® • , can be found by setting the first 
derivative of (3-80) to zero, vhich yields the following equation for 



R»« 






<XKl 






Ue'-B) J (3-81) 

For a sphere within an infinite plasma in i.iiich ionization is occuring, 
this is the radius at r4iich the effect of the spherical surface be- 
comes zero. 

The current density within this plasma is 



Q. r -77,r;D5CSC ] ■' 

xf4L_ - ■^^^<^frr{r-e)-j 

IfT^) 1-R'- J ' L J 



( 3 - 82 ) 



It is to be noted that it has two components. The first is directed 

toviard both surfaces from a radius within the plasma, and the second 
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is directed radially outward from the center, but is zero at both 
surfaces. 

For a spherical shell of radius R' with no internal sphere, one 
obtains 



with 




(3-83) 



J- JL JL 

Ds ■ R' (3-83a) 

This equation tos obtained from equation (3-80) by setting R = 0, and 
by evaluation of by LMlospital's rule 



C, ^iCo:i(^') ^ T|[ 

r->o p — ' ^ Ri ~ 



y}i 



(3-83b) 



where n^ is the density at the center. Figure 3-6 is a graph of this 
function. 



5. Prolate Spheroidal Geometries (the wire or probe) 

The influence of a prolate ellipsoid upon a plasna is now 
considered (see Appendix I for transf onnat ion equations and figure). 

In its limit, vihen the minor axis is small and the major axis approaches 
the length d, the effects of such a body can be used as ar approxi- 
mation of the effects of a finite wire or narrow probe of length d. 

Flamraer shows chat the i.»ave equation 

k^) W = 0^ (3-81i) 

equation (3-7) j expands in prolate spheroidal coordinates toj 




(3-85) 
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figure 3-^ Ro.c\i<a-l Densti^ -for Holbco Sphere 



► 
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in which ve have set 



(3"85a) 

and that^ by the method of separation of variables, solutions to 
(3-85) may be obtained in the form of the Lame products 



(3-86) 

The functions Sjnn(c,"^ ) and Rr^n (c,^) satisfy the ordinary differential 

= O (3_85a) 



in vhich the separation constants and rtl are the same in both 

equationso It is a combination of such functions (of the first and 
second kinds) that is given in Equation (3-13) as the general solu- 
tion to (3“81 i) 



a. Prolate Ellipsoid in stationary infinite plasma 

Before proceeding further with the prolate vave func- 
tions, one special case^ analogous to that of the sphere in the 
infinite stationary plasma, should be treated. By imposing the 
boundary condition that at large ^ the density^ ^ of a stationary. 



non 



■ionizing plasma be independent of and "7^ (so that = 0 



and '"^rpu = 0)j we can reduce equation (3-85) to 



vhich integrates tos 



(3-87) 
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( 3 - 88 ) 



vhere n 

o 

surface 



is the uniform density at infinity, and the density at the 
is 2erOo The current density is then; 







(3-89) 



To investigate the properties of this solution, ve shall use 
the identity 






(3-90) 



hy which Equation (3-88) becomes; 
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at very large distances (i.e,, ^ » 1), In (^+1) - ln(^ -1) 

and equation (3-91) can be expanded to 



(3-91) 



7t ^ 





(3-92) 



which is identical to equation (3-75) for spherical geometry, a 
result i-tiich reflects the fact that the surface "looks like" a sphere 
at large distances. In the vicinity of the surface, i.e., for 
vre let 



(3-93) 

(with 6^ such that the mean free path, X , is much less than ^ ), 
and expand, obtaining 




(3-9iia) 
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(3-9hb) 
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^9hc) 



By neglecting second order terms in ^ and using the expansion 

J^("X+ /) if X + ^ ^ > 

equation (3“9iic) can be reduced to 



( 3 - 95 ) 
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— 5e* * ^ 



( 3 - 96 ) 



Thus, for the density at the surface, ^ - 0 and therefore n = 0. Equa» 
tion ( 3 - 96 ) shovrs that the density is a linear function of the 
distance from the surface for short distances. 

At large distances, the current density can be seen from Equa- 
tion (3“92) to be an inverse sqiaare function of distance (as it is 
for a sphere); 






r 

At short distances, we obtain from (3-96) 



( 3 - 97 ) 



p. 



(3-98) 



vihich is a minimum at the center of the surface ( 7 ^ = 0) and approaches 

infinity at the ends (>| ^ t 1 ), This occurs because the area at the 

ends is zero, and thus apj' current flovring tot>>ards then constitutes 
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an * inf inite‘ current density. 

The current density at the surface is found fron (3-98) to be 



p ^ 

^^^hich again becomes infinite at the ends and is a ninimun at the 

center. Ve can obtain an estimate of the heating at the surface 

per unit length if we approximate the ellipsoid vrith a cylinder of 
L 

radius r “ ( - l)*‘and approximate the current density vjith that at 



the center = 0), 

£i 



r.StAA.jdLc^ 

/u 









4'n Aaic ) 

fo 



( 3 - 100 ) 



To find the behavior at siuall radii, vz let | = 1 obtaining; 

^ r-TFr-r— o-wi) 



^ (•/,”«.) [War f,.-) -^fcO 



which approaches zero as approaches zero, 

b. Prolate Spheroid in steady state plasma - ionization 
occuring 

Given an infinite plasma of homogeneous density ng at 
large distances and a prolate spheroid vihose surface is the surface 
, one can simplify equation (3-l8) considerably by consideration 
of boundary conditions. First, the requirement that the vtave function 
be finite at Tj s + i confines the solution to angle functions of 
the first kind since the angle functions of the second kind are singular 
at these points. If ve look for those solutions for vhich the den- 
sity^ n^is uniform at large distances for all f , then N is a function of 
^ and ^ alone, and m ® 0, Further, the prolate angle functions 
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reduce to the associated Legeadre functions when c vanishes; i.e., 
as the spheroidal insulator becomes spherical in shape. Therefore, 
vre can determine the subscript V in Equation (3-l3) by studying the 
Legendre function^, P^( ) because our solution must satisfy the limit- 

ing case of a hernispherical body. Only Pq( ) is positive for all Tj 
in the range -1'^ /j 1. Because n jnust be positive for all /y , the 
subscript must be 0. Consequently, v>e obtain for the density distri- 
bution 



n = K Sl? (c,^) *- br»‘ ( c , -5) 



( 3 - 102 ) 



where 



r ^ ^ 



(3-102a) 



(3-102b) 



If we employ boundary condition (2-99) ( >1 v^ll ■ 

ft 

ft 

It must be remarked that this solution does not satisfy the require- 
ment that the density be independent of at large distances. Hov;ever, 
for small values of c, and therefore for either very small or for nearly 
spherical bodies the dependence is negligible. The deviation is 
only for c - 0.5> but increases to practically 95^ for c = 5.0, 
c. Prolate spheroid and infinite vail 

Let us consider a semi-infinite plasma of uniform density 
n^ at large distances^ bounded only by a single infinite plane v-hich, in 
prolate spheroidal coordinates, is = 0. Further, consider an 
insulator (probe) V’iiose surface is defined by projecting into the 

plasma as indicated in Figure (3-7). The solution to the continuity 
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figure Prala^c Spheroid o.mA \r\^Yydc L)oJl 
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equation is developed in the same manner as in the previous section, 
except that in this instance the appropriate Legendre function is ^ ) 
because it has the property of being positive for all values of 
in the range 0 to +l and further of being zero at ^=0. Thus the 
density is 



77 = % So! L^o! ('s 0 ^ ^ ^ , 

where, again, 



%all “ 



_ M _ d 



Z4Vs. 



Ko’. 



B = - 



(3-103) 



(3-103a) 



(3-103b) 



The field and potential equations along the axis of the probe 
(i.e,, for >1® 1) arej 



(3-lOli) 



(3-105) 



fp ) - ^ Vf-b. 

f ;f- BKS'a,\>3 

vhere the prime indicates the first derivative, and 

- 4- £ --'0- g.. f + BRLVg,T^V 

^ (.Ro'’.(c,54 + 

6, Oblate Spheroidal Geonatries (circular lamina) 

In general, the solutions to the continuity equation in 
oblate spheroidal coordinates are the same as for the prolate with 
-ic replacing c and +i^ replacing ^ in the equations. Consequently 
these vtHI not be further discussed, except for the problem of the 

circular lamina in an inf init^ stationary, non-ionizing plasma. 
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a. Lamina in stationary plasma 



To stuc^ the influence in an infinite plasma of a 
circular lamina, the effect of an oblate ellipsoid is first discussed 
and then the limit as the semi-minor axis approaches zero. 

Proceeding as for the prolate case, by imposing the condition that 
at large ^ the plasma density be independent of and yj , and by sub- 
stituting -fi ^ for V vre obtain from equation (3-85) 




Integrating tx-.rice, ve obtain 



( 3 - 106 ) 



Letting 






?! =: w vp . 

71 ^ ) vre get 



7l^ r 

Using as the boundary condition 



hSll c^■^ 



( 3 - 107 ) 



(3-108) 



(W^ 



-O =■ C I ' Cjo) -tC^ 



subtracting this from equation (3-108) yields 

/> . -Jk— 

/.57l-fa^'(^o') r 



( 3 - 109 ) 



( 3 - 110 ) 



r - - 

’’ ^ y.57 1 -•WL'(^,') 

so that, for an oblate ellipsoid in an infinite plasma, 

\ h‘571 ■- 
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( 3 - 111 ) 



( 3 - 112 ) 



Vlien the semi=minor axis is zero, - 0, and the figure is a 
circular lamina of radius r^ ® c. Equation (3-125) becomes 



71 ’- 



7i^h,Z'(l) 



1.571 (3-113) 

The zone of influence of this amina, within vhich n ^r^n^^is given 



= <^* 725 ; 



(3-1 lit) 



i.e 



• > 



^ = (3-Ula) 

which is an oblate ellipsoid with semi-major axes r = v/^Tq and 
saai-minor axis z ^ r^o 

The current density in this plasma ( with c ) will bej 



Py - fe* 



This current density is a maximum at the surface of the lamina. 



(3-115) 



decreases as to zero towards infinity, and because the area of 

the "side” of tie lamina is zero, is infinite at the sides ( =0) of 

the “surface”. If we assume the current density at the center of the 
surface to be an approximation to the surface current density across 
the entire surface, then the recombination energy liberated per unit 
surface vull be approximately 






7TUi>U.ro 



(eV; •!§-■) 



I.S7I V- - ‘ ■ e 
xjhich, as for the sphere, is proportional to the radius. 



(3-116) 
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APPROXIMATE SOLUTIONS FOR INSULATED PROBES 

VTien an insulated probe is inserted into a dense plasma — that 
is, a plasma in which the mean free path is equal to or less than 
the probe dimensions— Langmuir' s probe theory £83 may not be appli- 
cable. Particle losses at the probe surfaces may cause appreciable 
effects on the particle density distribution, and hence on the probe 
potential. Under these conditions, of equation (2-71) can be 
appreciably large and consequently so may the effect of the in- 
sulator on the steady-state distribution. The approximation describ- 
ed in Section III-A-2 has been employed to calculate this effect 
for a floating probe in a partially ionized gas. Two different 
geometries are considered. The first is a prolate llipsoid 
inserted radially into a cylindrical discharge, with the sensing 
element on its tip. The other is a cylindrical probe aligned parallel 
to the discharge axis. For the latter, ve are able to account for 
the effects of a longitudinal magnetic field. 
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A. The Prolate Probe 



We do not have an exact solution of the continuity equation 
(2-55) for a cylindrical or ellipsoidal probe projecting from a 
cylinder ijall into a plasmao However, we have seen that there 
is a solution for the density distribution in a long, narrow 
plasma (equation (3“66), In the absence of a magnetic field, it is 



n = rio ^o(J^ f), 



ih-i) 



in which 



IT Z40S - (h-la) 

" R . 
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Further, we have obtained a solution for a prolate spheroid attached 
to an infinite wall (equation (3-102)) 

n = nJ So, 0 B (u-s) 

Based on the assumptions outlined in Section III-A-2, we shall 
approximate the effect of a narrow cylindrical or prolate ellipsoidal 
probe by joining the tw models described above as indicated in 
Figure ii-l. The infinite vrall is tangent to the cylinder, and the 
2-axis of the probe ( >^ = 1) lies along the diameter of the cylinder. 
We assume that the presence of the probe does not significantly alter 
the ioni 2 ation or diffusion rates, so that equation (3-103a) becomes 



d IZ 



0i-2a) 



vrhere i js/ , the length of the probe (if the minor axis is much 
shorter than the major axis) and R is the radius of the cylinder. 

No effort will be made to describe the density throughout the 
cylinder. We shall merely attempt to predict the density distri- 
bution along the ® 1 axis of the probe (i.e. , along the diameter 
of the cylinder). It is this density distribution v'hich is res- 
ponsible for the potential distribution sampled a sensing element 
at the tip of the probe. In accordance with the procedure outlined 
in Section III-A-2, we readily obtain the following condition for 
determining the position along this axis of the dividing surface be- 
tvreen the two models 



1%1 



:l^l 

> bylii 



•Probe ’^'Cylinder (U-3) 

Equation (U-3) was solved graphically because sufficient published 
data on the derivatives of the spheroidal wave functions was not 



available. 



IS 



Because of this paucity of information on the spheroidal vave 
functions, it vras also necessary to limit the study of this problem 
to the limiting case where X is so small that the boundary con- 
dition (2-99) (n 1 , ® 0) i^ay be employed. For this problem, B of 

Ta’cI 1 1 

equation (li-2) is readily obtained! 




The normalization constant, nj, of equation (U-2) vras determined 
in relation to n^ of equation (h-1) by the graphical solution of 
equation (U-3). 

A study of the solutions, and of the physical properties of our 
model, reveals that this model can be justified as an approximation 
to the particle distribution in the vicinity of the probe for probe 
penetrations varying from zero to about half the cylinder diameter. 

For short probe lengths, particle losses to the adjacent regions of 
the infinite vail (vrhich, for the short probe, is a reasonable approxi- 
mation of the adjacent cylinder vail) is the predominant factor in 
establishing the particle distribution in the vicinity of the probe. 

For instance, for a probe length of 0.2R, the parameter c of equations 
(U-2) is 0.U8. From Figure II-l (Appendix II) v'e see that the pro- 
late angle function is very nearly that for c = 0, vjhich 

corresponds to zero probe penetration. Consequently, the effect of 
the short probe is very slight in comparison to the effect of the 
vail, as vre expect from consideration of the surface which vrould be 
the cause of most of the particle losses from the region around the 
probe. However, for deeper penetration, the angle function changes 

appreciably, indicating that the surface of the probe exerts a more 
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dominant influence. Again, this is expected since particles in the 
vicinity of the tip of the probe are more likely to diffuse to some 
portion of the probe than to succeed in reaching the vail. Thus, 
even though the infinite vrall is no longer a good approximation of 
the cylinder viall for the deeper penetration, the increased probe 
surface losses justify using the resulting particle density distribu- 
tion functions to describe the plasma near the end of the probe. 

TfJhen the probe penetrates more than half-way through the plasma, 
the influence of the opposite walls of the cylinder on the particle 
distribution near the end of the probe become increasingly more important, 
and our model breaks dom. 

Conversion from prolate coordinates to probe dimensions is accomplished 
vrith the following equations. 



probe lengths 







probe radius ( 



semi-minor axis): 








is the equation of the probe surface. 



Flammer [63 t'jas the source of the information on the spheroidal 
irave functions. Data not calculated or graphed in this reference is 
presented in Appendix II. 

The graphs on the following pages were computed using the proce- 
dures outlined in this section and subsidiary equations developed in 
preceeding sections. 

Figure U-2 is a graph of density versus distance for probes of 
several radii but equal depths of penetration. Note that the magnitude 
of the dependence of the probe influence on the probe radius is 
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relatively small, an order of magnitude change in probe radius 
causing an increase in the region of influence of only ^0^. 

Figure li*»3 is a graph of the density distribution for various 
depths of penetration of a probe of radius rp = 0.02R (vhich, for a 
2.5 cm discharge, corresponds to a 0.5 mm probe). 

Figure li-ii shows “effective probe length" vs. true probe length 
for the 0.02R probe. The effective length is based on the inter- 
section of the two models, and is approximately the point of the ori- 
ginal particle distribution that the probe samples. The curve marked 
X 2r is the effective probe length when the probe does not 

disturb the plasma, and corresponds to the true probe length. The 
broken segment of the curve is an estimate for that region in vhich 
this model breaks down, as outlined above. 

Figure li-5 is a grajii of the potential disturbance associated 
with various depths of probe penetration. This potential disturbance 
is defined as the potential drop betvreen the point the plasma samples 
(the effective probe length) and the sheath edge. A sheath depth of 
O.OIR was assumed for this calculation. These potentials can be 
used with the effective length data to analyze experimental results 
or to predict them from the theory. 

Figure li-6 represents such a potential prediction. It is the 
prediction of the potential measured by such a probe assuming that it 
actually measures a difference between the potential at the maximum 
density legion of the plasma (zero potential) and the potential at 
the sheath edge at the tip of the probe. 

A very significant prediction of the above theory is that the ise 
of an insulated probe in a plasma of diort mean free path will cause 
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a definite asyranetry even in an originally symmetric plasma, 
these results are not applicable directly to a plasma in the 
of a magnetic field, some of this effect is still expected. 



Although 

presence 
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B. Qylindrical Probe Parallel to Axis 

Figure h~l is a top view of a possible insulated probe con- 
figuration for measurements in the presence of a magnetic field. The 
cylindrical probe (enlarged in the figure) is supported at its ends 
and long enough that the effect of the end supports are negligible 
in the region at the center of the probe. The sensing element (or 
elements) are located at the center of the probe on the common 
diameter of the probe and plasma. 

There is only one alignment of the probe for which an exact 
solutionof the continuity equation (2-57) is possible: when probe 

and cylinder are co-axial. As previously developed (equation 
(3-72)), the solution is 






r < n I? 



Oi-7) 



with r the probe radius, R the cylinder radius, and n* the 
normalization constant. If A is not very much less than the probe 
dimensions, boundary condition (2-71) should be used at the probe 
surface. However, equation (2-99) (*^^^2,1 * often still appli- 

cable at the cylinder v^all, and was used for our calculations. 

Using equation (2-71), which is 



11 


luatC. 




(ii-8) 


sheath yields, from (li-1) 
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(li-io) 



from which we obtain 




Next, employing equation ( 2 ^ 99 ) at the outer cylinder, we obtain 



= CNX&Ji) 

Equations (U=10) and (u“li) are sufficient to determine both C and 

and they are most easily solved by successive approximation, 

‘ • - , 

Table U“1 lists the results for these parameters for various K between 
zero and one, using the following probe dimensions and mean free path; 



r ® 0.02R 
P 

A = 0.02R , (li“12) 

(again this corresponds to a 0,5 mm probe in a 2.5 cm cylinder). 

Even when there is no recombination at the probe ( (X, =0), the rate 
of ionization for the steady state must be greater than when the probe 
is not present because the probe still occupies a finite \o bone in 
the cylinder, 

When the ,probe is off=axis, it is necessary to employ the 
approximation described in Section III“A“2. Equation (ii“3) is still 
applicable for determining the separation of the zones along the 
common diameter (the only location in >hich we inll be interested). 

In order to evaluate vicinity of ,the probe for various 

probe positions, v»e interpolate between the calculated value when the 
probe is coaxial and the value it vtould have when the probe is at 
the cylinder i-?all (and has negligible effect on the eigenvalue because 
of the low density there). This latter value we assume to be th«~ 
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familiar 2ohO$/R<‘ Assuming that this parameter would depend on the 
density of the plasma in the position in t^hich it is placed, we 
obtain the following relationship 




(li-13) 



where 

From the values thus calculated, a mean was calculated for each 
(see Table ii=l), and this value ■nas used in the computations of the 
density distribution in the vicinity of the probe. Tlie use of a 
more precise value is not justified by the nature of the approximation. 
For the density distribution in the regions controlled by the cylinder 
vfalls, the value of J for the cylinder alone vras used ( 2 ,Ii 05 /t?) 
because the exact effect of the probe on this distribution could not 
be accurately estimated. The degree of error resulting from these 
approximations can be inferred from a comparison of the exact solution 
for the coaxial probe and the approximate solution shovn in Figure 
U“>lii. As expected from the approximation model, the maximum error is 
in the region where the two models are joined (where the positional 
error is O.IR and the magnitude of the approximation is 12^ greater 
than that of the exact solution). The magnitude of the error de- 
creases towards either surface, but the percent error remains approxi- 
mately constant. Consequently, the potential determinations (vhich 
are based on a logarithm of a ratio of densities) vrill be reasonable 
approx ima t ions , 
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Several distribution curves have been obtained by using the 
procedure outlined above? (1) Figure ii~8 is a graph of the exact 
solution for a 0o02R® co=axial cylinder for various oC ; (2)Figure 

is a graph of the exact solution for co-axial cylinders of various 
radii at a given (3) Figure li-lO is the density distribution for 
various (ivi for a 0.02R probe in an off-axis positionj and (h) Figure 
li-ll is for a given c< and various probe positions for a 0»02R® probe. 
Effective Probe Position vs, actual probe position is indicated 
in Figure ii-12 and potential correction to be applied to probe readings 
in Figure li-13s we nave assumed the sensing point to be on the side 
facing the center cf the plasma when the actual position is less 
than |-R from the wall. Potential curves are given in Figures h-lh 
and h-l5» The former one is for one contact, oriented as indicated 
above. The latter is for tw*o contacts, one on either side of the 
probe along the common diameter. Figure ij-13, the potential correction 
curve, is just the potential difference between the sheath edge and the 
point of intersection of the zones of influence for various p<i and 
probe positions. 

Figure li-l6 is a curve of heat liberated at the probe for 
various values of probe radius when the probe is aligned with the 
center of the cylinder^ and ii-=17 is for various values of o<: , It is 
assumed for these calculations that the fraction of particles which 
diffuse to the sheath reach the surface and recombine, liberating 
(eV^ -r ) of energy. 

The results of the above theory indicate that the cylindrical 
probe should have considerably less effect on the particle distribu- 
tion functions than the prolate probe in a given cross-section of 

the plasma. This wrauld be expected from the decreased surface 
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involved for the cylindrical geometryo Consequently, there is less 
assymetry produced by the probe, and the potential curves, particularly 
vith contacts on both sides of the cylindrical probe, should be norc 
synmetric. Hovever, the effect of neasuring a potential at an 
"effective distance" causes some flattening of the predicted curve, 
Figjre h~9 indicates that the density distribution in the vicinity 
of the probe is relatively insensitive to the probe dimensions. 
Consequently, the potential drop associated vtith this probe configura- 
tion is relatively constant for various probe radii. 
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TABLE li-l CYLINDRICAL PROBE PARAMETERS— EXACT VALUES OF 



values vs® oC. 



V'S. oL for probe at center, and mean 
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Mean Values 
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Co Evaluation of Probe Keasurer.ents 



In a dense plasna, in vhich the mean free path is of the 
order of magnitude of or less than tlie probe dimensionSs strong 
deviations from Langmuir's theory can be expected. These deviations 
depend on the size and character of the probe. By using the 
approximation described in this thesis, one can calculate these effects 
in terns of an "effective probe length" and a ‘’potential disturbance." 
These tw3 parameters can then be used to analyze or predict the 
experimental results in the mariner just employed for the prolate and 
cylindrical probes. 
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MACiWSTIC FIELD ETID EFFECTS 

Let us nov investigate the effects on the pax'ticle density 
distribution functions of the decrease in magnetic field intensity 
at the end of a coil vith attendant intersection of field lines vith 
the container v,»alls. The effect expected is that diffusion along the 
field lines, being greater than that perpendicular to the I’alls (radial 
diffusion), vould cause the density distribution near the ends of the 
field coil to decrease more with increasing radial distances from the 
axis than predicted from radial diffusion alone* Particle losses by 
longitudinal diffusion to the end of the field coil, lAere sc<me of 
the field lines intersect the container walls, would be responsible 
for this effect, which v’ould manifest itself physically as a con- 
striction near the ends of the field coil. 

An exact solution of the continuity equation for this configuration 
(Figure 5“l*"note the enlarged radial scale) could not be found, and 
consequently the approximation described in Section III“A»2 was used. 
The plasma is divided into zones in vhich the diffusion ^■ould be 
primarily radial or primarily longitudinal (along field lines), 

ITie latter zone, in the vicinity of the intersection of the field 
lines with the plasma container, is a zone in which the particle dis- 
tribution is controlled by the loss of particles diffusing along the 
field lines and recombining at the region of the wall where the lines 
intersect the wrall. The size of the region of parallel diffusion 
increases, of course, wdth increasing magnetic field. If wre neglect 
the curv'ature of the fielc lines and use the boundary condition 



102 



ve obtain for the density distribution along a field line 



n(z) = n 







(5=1) 



vfherc n* is a normalizing factor to be determined j the z=axis is 
the cylinder and magnetic field axis, and the origin is the inter~ 
section of the field line and the v,'all» TTiis solution cannot be 
rigorously derived, but is justified as an approximation consistent 
VTith the accuracy of this calculation. Comparison vith the solution 
for a bounded cylinder in a homogeneous field, equation (3=62), shovs 
the same z-depcndence for the density distribution, 

A second zone in this plasma is that in the vicinity of the center 
of the field coils in vhich vfall losses due to radial diffusion to 
the rialls control the density distribution. Here 



where is the density in the center of the cylinder. 

Then, as a first approximation, vre shall assume that the ioniza= 



significantly from that of the corresponding stable radial configura= 
tion. In higher approximations we would have to calculate the change 
in this quantity due to induced changes ir electron temperature by 
the different fields, but v^e shall assume these to be rcgligible. 
Consequently, 




tion coefficient, , in the zone cf parallel diffusion wilt not ciffer 



S>' ^ ]) 

Therefore, we obtain from equations (5“l) end (5=2) 




( 5=1 



103 



The surface separating the two models (equation (3*hl)) is that at 
which 




5 - CC^(^ -) 



P, 

Xi^r) 






Substitution from equation (5“lj) reduces this to 




Figure (5-“?) is a nonogram from which the z-coordinate of the 
separation surface, z, can be obtained as a fraction of z^, using 
r (as a fraction of cylinder radius R) and 




as arguments. From equation (^=7), it is apparent that **a'* decreases 
(with B) along the cylinder axis (indicated on Figure (5“1)) beyond 
the end of the field coils. Consequently, vhen using the nomogram 
to determine Zq, it is necessary to use the value of "a" at the 
junction of the tv^o models. This requires a trial=“and=error procedure. 
In order to illustrate this method of approximation, and to 
estimate these effects, a systan writh the following dimensions v»as 
selected; 

Plasma container radius; 2,? cm, 
magnetic coils; wridth ~ 1 $ cm, 

spacing «= 5 cm, 
number “ 10 

inside diameter - 2 $ cm, 

outside diameter - 50 cm, 

overall length - 195 cm, 
lOii 



The axial field of the above array of field coils vas computed with 
a digital computer (figure 5-3) > using the thick coil equations from 
reference [20} » The author states that they agree to within a 
few' percent of measured values, and for the region within the coil 
array, experimental curves from a very similar array were available 
and agreed very wall* The field was assumed constant through any 
cross-section of the cylinder, which is a good approximation because 
of the small radius of the cylinder. Figure 5-1 shows the field lines 
as computed from the field strength data. 

Figure (5-1) also shows the separation of the tw zones of 



diffusion for three different orders of magnitude of namely, 

/O 






_ tOOO-l 
U 






V-} - OOP - ' 

0.-100 

(5-8) 

It is to be noted that for B (a ® 10) the effect of the 

/V'' 

parallel diffusion is very small, and is centered at the end of 

the coil array. For « 31.6) the effect is appreciable, 

and extends more than half-way into the cylinder. Figure (5-ii) 

is a density profile of the resulting particle distribution; Figure 

(5-5), selected cross-sections from Figure (5-^4)! and, Figure (5-6) 

is a plot of constant density curves. 

One feature of the density distribution which must be discussed 

is the question of detemining the density distribution of the central 

zone v?hen it does not extend to the trails but rather terminates at 

the edge of a longitudinal zone. A longitudinal current (between 
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electrodes at either end of the cylinder) is responsible for the 
ionization uithin the plasna, and this current must be constant 
through the plasma. Further, we have assumed that the ionization 
coefficient, p , is constant. From this it follows that the electron 
temperature is constant, and therefore the longitudinal electric field 
must be approximately constant. But if the field and current are 
approximately constant, then the total number of particles in a 
cross-section must be constant. Further, the plasma in the radial- 
diffusion region vhich is not controlled by the cylinder boundaries 
because of the intervening region of parallel diffusion is exj^ected 
to maintain its Bessel distribution (thus satisfying the diffusion 
egruation and not requiring a different ionization coefficient). As 
a consequence of the constant field and current restrictions, however, 
the total cross-section particle density (particles per unit length)^ 



is a constant, equal to the total number of particles per unit length 
at the center of the field coil. Therefore the normalization of the 
Bessel distribution for the center of the plasma in regions i.here the 
distribution at greater radii is controlled by longitudinal diffusion 
must be so determined that equation (5-9) is satisfied. 

One effect vhich ve ignore in the previous formulation is that 
due to the small radial component of the magnetic field at the end of the 
coil. This component is strongest at greater radii, and decreases 
the mobility of the particles In the longitudinal direction. Con- 
sequently, it vculd contribute further to the constriction of the 




( 5 “ 9 ) 
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plasma at this point, V'e neglect this influence, assuming the 
curvature of the field lines to be very anall (recall that the radial 
scale of Figure (5-1) is greatly enlarged). 

It should be noted that there is a discontinuity in the direction 
of V n along all the dividing surfaces betv?een zones, and a discon= 
tinuity in particle density, n, along the a » 31,6 boundary from 90 cm 
to 130 cm and the a ® 100 boundary from the center to 161^ cm. The 
former discontinuity is a consequence of our approximation. The 
density discontinuity occurs because the dividing surface lies along 
one of the field lines, and is a result of the fact that the sine dis= 
tribution for the longitudinal model cannot be made to match the 
normalized radial distribution along such a line except at the end 
which is nearer the end of the coil. The transition regions are 
centered along these dividing surfaces, and are the region of greatest 
uncertainty. Along those surfaces at vihich the particle density, n, 
is discontinuous, this approximation is even more uncertain. Because 
of this, the a = 100 description, with a discontinuity in both 
and n along its entire length, must be much too crude to be valid, 
and further attempts to fit this approximation to the effects of 
such large fields are not attempted. The description for a = 31,6 
is believed to be a valid '‘order-of-magnitude" approximation. 

In Figure (5-7) the heat liberated at the walls ty the ions and 

electrons carried to them is indicated as a function of position along 

the cylinder. There are three curves; all calculated for a = 31o6s 

(1) the v’a.ll heat loss for normal ambipolar diffusion without a 

magnetic field, (2) the -rail heat loss for radial diffusion with a 

magnetic field, accounting for the increase in radial diffusion t«hich 

accompanies the decrease in magnetic field outside the coils, and 
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(3) the estimated heat loss for the magnetic field, taking into account 
the longitudinal as well as radial diffusion. This latter curve is a 
very crude approximation based on the results predicted by this model 
with adjustments to compensate for the errors introduced the 
discontinuities at the separation surfaces. The figure shoTf»s that 
the heating at the irfall vrithin the region of parallel diffusion and 
within the coil array is greatly decreased because the particles are 
transported to the end of the cylinder instead of directly to the Tr’alls. 
However, outside the coil, v»here particles from within the coil are 
carried by the longitudinal diffusion, the effects of such additional 
heating is masked by the heating due to the increased raidal current 
which flows because of the lov^er magnetic field intensity. 

The values of magnetic field B for v.hich this effect can be 
approximately described by this model are those v^iose order of nagni=> 
tude corresponds to a <»31.6. Smaller fields (a ^£:10) will not cause 
an appreciable end-effect constriction. Larger fields ( 0 - 2^ 100) 
are not describable by this model. For a plasma which fits the sample 
data given in section 3-B-l, this is seen from equation (5-8) to be 
approximately 



It must be recalled, however, that the above description is based 
upon the vails being perfect insulators. Electron diffusion, pre- 
dominately parallel to the field lines because of the small electron 
Larmor radius, is restrained ly the build-up of a negative v-all charge 
in the vicinity of the magnetic field line intersections with the 
cylinder vralls. Similarly, the ion diffusion, which would be predominantly 
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radial because of the larger Larmor radii, would be opposed by a large 
positive wall cliarge within the coil. In this manner, the conditions 
of equal electron and ion radial current and equal longitudinal current 
are net. However, it is known that insulators have a surface conduc- 
tivity vtiich depends on the material and the temperature. The effect 
of this conductivity is to permit an exchange of carriers along 
the rail. Consequently, we obtain diffusion vhich is sonevtiere betvreen 
Allis diffusion [1^3 ai^d our quasi-ambipolar diffusion. The vtall 
charge build-up is avoided by the exchange of ions and electrons 
between negatively and positively charged zones. This exchange results 
in recombination along the length of the cylinder and would probably 
invalidate the conclusions of Figure (5-7) if the currents were 
appreciable. Further, if this surface conductivity were sufficiently 
high, the diffusion would approach Allis diffusion, under vhich condi- 
tion the ratio of the parallel and perpendicular components of the 
diffusion tensor would no longer be 1 + but vrould be 

1 * B? Thus the field strength «diich would produce a con- 
striction describable by this model would be reduced to 



The description of the density distributions would not be affected 
by the change in diffusion process since the only effect on the 
diffusion is a change in the magnitudes of the components of the 
diffusion tensor. In practice the diffusion would be somewhere be- 
tween the extremes of no surface conductivity and perfect conductivity. 
Therefore, the value of B to vhich these calculations correspond vrould 
intermediate betvroen 0,li veher/n^ and h veher/m^f depending upon 




(5-11) 



the surface current conducted by the pla3na container walls, 
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It should be noted that this description predicts negligible end effects 
for the plasma described for fields below about 0.1 wetaer/m^. For 
particles with higher mobilities, however, this effect could occur in 
smaller magnetic fields. 
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APPEM)IX I - Curvilinear Coordinates; Transformation Equations 



Because of notation differences between references, the trans- 
formation equations for all coordinate ^sterns employed in this thesis 
are presented below; 



1, Cylindrical Coordinates {r,5^,2') - Figure I-l 
X = r cos^^ 

y = r sing( (l-l) 

z = z 

where r 2. 0, -oo<z< oo 



2. Spherical Coordinates (r,9,<^) - Figure 1-2 
X = r sin9 cos^ 
y - r sirQ cos0' 
z = r cos9 

where r>0 0S?fd2-n O^e^TT 



3. Prolate Spheroidal Coordinates (/) - Figure 1-3 



vrhere -1 ^ < 1 1~^'<cd 0 ^ 2tj 

In the prolate spheroidal system the surface ^ = constant >l is an 

elongated ellipsoid of revolution vuth major axis of length and 

■ 2 . J- 

ninor axis of length d ( J -1)2, The degenerate surface y ® 1 is the 
straight line along the z-axis from z = -^-d to z = 4;|-d, The surface 
1>/| = constant <1 is a hypei'boloid of revolution of two sheets with an 
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asymptotic cone vhose generating line passes through the origin and 
is inclined at the _^|3 “ cos“'‘-/| to the z-axis. The degenerate surface 
|7|| = 1 is that part of the z-axis for which |z| < ^d. The surface 

= constant is a plane through the z-axis forming the angle ^ with the 
x,z-plane. 



ii. Oblate Spheroidal Coordinates (>|. f .s') - Figure I-I4 

^ “ "I &' -Mil ^ 

t where -1£>^^1 Oi^^oo 0i?f£-2lT 

In the oblate spheroidal system the surface )'’^| ^ constant > 0 is a 

T im i- 

and minor axis of length d|^|. The surface =0 is a circular disk 
of radius a - ^d which lies in the x,y-plane and is centered at the 
origin. The surface H - constant ^ 1 is a hyperboloid of revolution 
of one sheet with an asymptotic cone vhose generating line passes through 
the origin and is inclined at the angle ^ ^ cos^^/j to the z-axis® 

The degenerate surface j>j|=l is the z-axis. The surface 7 ^ =0 is the 
x,y-plane except for the circular disk ^= 0 , The surface = constant 
is again the plane through the z-axis leaking the angle ^ with the x, 
z-plane. The transformation from prolate to oblate coordinates is 
made by replacing ^ by i’^ in the prolate equations. 



Reference for spheroidal systems; Flammer, Spheroidal Wave Functions 6 
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APPEI'DIX II 



SPHEROIDAL ViAVE FUWCTIOWS 



The particle density distribution in the vicinity of a prolate 
spheroid is describable by the prolate angle function of the first 



kind, and the prolate radial functions of the first and 

second kinds, Rci and > respectivelyo 

Flannier [6} has tabulated So(Yq?|) (See Fig. II-l, which is 
taken from his data) and the radial functions and j^oi 

for different values of c and the independent variables and , 
The range of ^ the radial functions in his tables is, however, 

quite restricted ( ^ =1,0C^, 1.020, l.OUi, and 1.077) and consequently 
it TOS necessary to compute further tables of the radial functions 
(Table II-l and II-2). 

To compute Roi 0^,^) and f\oi the following expansions 
in terms of spherical Bessel functions and Legendre Polynomials were 
obtained from Flamner's general equations, (li.l.l^) for (c^^)and 

(U.2.6) for 

1 



(C, ^ £'ar (Il-l) 

nd ' 

where from Flammer’s equation (U.1.8) and (U.2.5b)j 

-jr ( f) ~ 



Il-la 



II=-2a 



^0 1 Lw ) - 3> I 

and where the prime over the summation sign indicates that the 

summation is over only odd values of r, and where ('I (c) and d|,|',(c) 
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are e>pansion coefficients. Flamraer tabulates these coefficients 
to five significant figures. 

The prolate radial functions of the first kind, K’oi'CC,^) (Table 
II-l), were calculated using a desk calculator and the Mathematical 
Tables Project Tables of Spherical Bessel Functions Ql] . The figures 
have been spot checked, and are believed correct to plus or minus one 
in the last digit. 

The Prolate Radial Functions of the second kind, ^^^^(c^)(Table II-2) 
were calculated on the C02-l60li Computer using equation Most 

of the Legendre Functions required for the expansions were also genera- 
ted on the computer, using the following equations given by Jahnke 
and Emde m . 




h 3 S---(2n-l) 



7, 7ih-' ) v)-^ 




Because of the slow convergence of the equations (II -3) for large 
values of ^ only those functions for '^<1,2 were computed. The 
values for ^ > 1.2 were taken from Tables of Associated Legendre 
Functions £j.2] . The Values of Qq for J < 1.2 were calculated on 
a desk calculator using the following relation from Jahnke and Emde [ 103 . 



Qjf) = cdrYf ) 

and tables of logarithms. 






11=5 



The radial functions are tabulated here 

for values of c between 0.1 and 2.0 and for values of J between 1.0 
and3.6 and 1.0 and 2.0 respectively. These functions are graphed in 
Figures (II-2), (II=3)> and (Il-li). The tabulated values are correct 
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to the number of digits shown, except that all five digit 
results are considered correct to plus or minus one unit in 
the last digit since Flammer's coefficients are listed to only 
five significant digits. The decreasing number of significant 
digits for higher values of c and ^ is a consequence of the 
decreasing rate of convergence of the series (II-2) for rela- 
tively large c and J , The number of significant digits v»as 
determined by inspection of the terms and generous estimates of 
the maximum magnitude of the neglected terms. The author has 
endeavored to err on the conservative side in reporting these 
values, and believes all digits reported to be correct. 



126 



J 

I 



I 

I 

I 



•V 



TABI£ II- 1 

PROLATE RADIAL FUNCTIONS 
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